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SUBSTITUTE SPECIFICATION 



METHOD FOR PACKING AND FILTERING GEOPHYSICAL EVENTS TO 
RETRIEVE THEREFROM DATA ON THE TYPE OF SUBSOIL 



ABSTRACT 

The invention relates to a method for packing and filtering geophysical events read 
on multi-domain records, with distribution of said events into families each having a 
particular geophysical significance : iso-offset or iso-incidence angle data cubes, elastic 
parameter cubes derived from a joint stratigraphic inversion and the like, in order to 
extract therefrom data on the type of subsoil, using a multivariate statistical technique. 
The method essentially comprises forming, by combination of geophysical variables, 
synthetic variables much fewer in number, which are obtained by constructing an 
orthogonal vector base in each of the sets of analysis constituted by the data of each of 
the families, wherefrom is derived the formation of an orthonormal vector base (new 
attributes) for filtering and describing said geophysical events. The invention is 
applicable to geological interpretation or to interpretation of an underground reservoir, 
from seismic measurements or logs for example. 
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METHOD FOR PACKING AND FILTERING GEOPHYSICAL EVENTS TO 
RETRIEVE THEREFROM DATA ON THE TYPE OF SUBSOIL 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 
Field of the Invention 

fOOOn The present invention relates to a method for packing and filtering 

geophysical data in order to retrieve therefrom information on the nature of the subsoil. 

This methodology is applied to packing and filtering sets of multi domain 
geophysical data, i.e. cither pre ■ stack seismic data from a 3D seismic survey grouped 
into a certain number of data cubes corresponding each, for example, to the same offset 
(seismic trace cubes referred to as iso offset cubes) or to the same angle of incidence 
(referred to as iso incidence angle cubes), or seismic attribute cubes obtained by seismic 
amplitude processing such as, for example, elastic properties cubes from a joint 
stratigraphic inv e rsion of pre stack seismic data, these cubes being obtained by 
implementing the method described for example in patent FR 2,800/173 (BP 1 ,096,270) 
filed by th e applic ant. It then consists in shap ing the pr e- stack seismic information 
contained in seismic trace portions formed by exploration of an underground zone, by 
means of new representative attributes which arc thereafter used for geologic 
interpretation or interpretation of an underground reservoir. It can also apply to packing 
and filtering of pre stack or post stack 3D seismic data from repeated surveys (4D 
seismic method) in order to interpret the fluid content variations of a petroleum 
reservoir under production. This method is also applicable to packing of log data of 
several physical natures by forming, from initial m e asurements, new log attributes in 
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order to improve interpretation of the log data in terms of lithologi c an d p ctro physi cal 
characteristics of the reservoir. 

BACKGROUND OF THE INVENTION 
Description of the Prior Art 

[00021 Geophysical measurements are conventionally used to provide 

complementary information, in relation to direct observations in wells, on the geologic 
properties variations of underground formations. Seismic interpretation is for example 
based on the analysis of attributes in order to retrieve these geologic properties from the 
seismic information, in a zone of interest of an underground formation (reservoir zone 
for example). Another example is log interpretation where the various measurements are 
simultaneously analysed analyzed to deduce therefrom the petrophysical characteristics 
of the underground formation studied. 

f0003] A standard approach for seismic or log attributes analysis involves the 

use of statistical pattern recognition and estimation methods in order to calibrate the 
pertinent attributes with the available information obtained in various wellbores. 

TO0041 F or example, paten t FR French Patent A - 2,768,818 and corresponding 

fUS[- ] Patent 6,05 1,651} filed by the applicant assignee describes a method allowing te 
detec tion of the potential classes in a population of events related to the physical 
properties of a complex medium such as the subsoil, located from data obtained by 
subsoil exploration or by in-situ measurements (events of geologic, geophysical nature, 
events linked with fluid production data, etc.), these events being located by points in a 
multi-variable space and defined by an attributes vector. 
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rOOOSl Another example is given by patent FR - French Patent Application 

01/05,675 filed by the applican t assignee , which describes a method for facilitating 
monitoring with- over time of the evolution of an underground zone by compared 
analysis of a certain number n of seismic records successively obtained after n 
successive 3D seismic surveys (4D seismic method), wherein a pattern recognition 
technique applied to the whole of the seismic events of several surveys, considered and 
analysed analyzed simultaneously, is used. 

[0006] Patent EP EP Patent A-671,017 and corresponding f US[- ] Patent 

5,638,269) filed by the applicant assignee can also be mentioned, which describes a 
method allowing to-establish ing of a connection between geologic data obtained by 
coring or logging in wells and seismic data obtained by means of a seismic exploration 
survey, which is based on a statistical calibration technique with statistical calibration 
obtained by bringing together the local geologic data (measured in wells) and seismic 
attributes read on seismic traces obtained in the immediate vicinity of each well. 

1*00071 The seismic attributes that are subjected to such an interpretation are 

calculated in the post-stack or in the pre-stack domain. 

[0008] In the post-stack domain, the conventional attributes are calculated from 

the amplitudes at the reservoir level, or from the impedance P estimated by means of a 
stratigraphic inversion. 

[0009] In the pre-stack domain, the number of pertinent attributes can increase 

considerably. In this case, the attributes are either the amplitudes at different offsets or 
incidence angles, or parameters obtained from a pre-stack joint stratigraphic inversion 
(impedances P and S, density, product of the Lame parameters with the density, etc.). 
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fOOIQI Concerning the interpretation of logs, the situation is similar^ jo a very 

large number of measurements being available, hence a large number of attributes 
allowing to characterize characterization of the logs, these attributes being often 
organized in families of equal physical significance (for example family of the 
resistivity attributes, family of the attributes associated with the natural or induced 
radioactivity of the formations, etc.). 

fOQ11] A method referred to as principal components analysis (PCA) method, 

well known to the man skilled in the art, can be applied to these attributes in order to 
analyse analvze the relations existing between them and to reduce the number of 
significant attributes. The principal components extracted from the PCA are the new 
attributes : they define an orthogonal or non-orthogonal (rotation) vector base and are 
linear combinations of initial variables. They can be used in statistical pattern 
recognition algorithms and correspond to a multivariate filtering of the initial seismic or 
log information. Various applications of the PCA method are for example described in : 

- Dumay, J., Fournier, R, 1988, "Multivariate statistical analyses applied Statistical 
Analyses Applied to seismic facies recognition Seismic F acies Recognition " 
Geophysics, 53, 1151-1159, 

- Hagen, D., C, 1982, "The application Application of principal components analysis 
Principal Component Analysis to seismic data sets Seismic Data Sets ", Geoexpl., 20, 
93-111, or 

pat e nt PR French Patent A -2,772,138 and corresponding (USj-] Patent 6,345,108) 
filed by the applican t assignee . 

[QQ12] However, PCA does not take into the account of the notion of attributes 

group and, consequently, it is often difficult to give a physical interpretation of the 
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principal components, or to relate them clearly to the initial attributes, especially if they 
are organized in groups of equal physical significance. 

fOQ13] The analysis referred to as canonical analysis is also a statistical method 

well known to the man skilled in the art, which allows to-establishing the relations that 
may exist between two sets of variables in order to know if these two sets describe the 
same properties. This method is for example described by : 

- Hotelling, H., 1936, "Relations b e tw ee n two sets Between Two Sets o f variables 
Variables ", Biometrika, 28, 321-377. 

[•QQ141 An example of application of canonical analysis is for example described 

by: 

- Fournier, F., and Derain, J.F., 1995, "A statistical methodology Statistical 
Methodology for d eriving reserve* prnpwtipr Deriving Reservoir Properties from 
seismic dat a Seismic Data ", Geophysics, 60, 1437-1450. 

r00151 This method is limited to the study of two sets of variables and does 

therefore generally not apply to the multicube seismic information or to the 
multidomain log information. Furthermore, it defines synthetic variables in each 
subspace associated with the two sets, and not a single vector base allowing to-packing 
the whole of the initial variables and to describe the relation between the two sets. 

f 001 61 There are several methods allowing to gonoralizo generalization of the 

analysis referred to as canonical analysis b y extending the analysis to more than two 
sets of variables. Different aspects of the prior art in the field considered are described, 
for example, in the following publications : 
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- Horst, P., 1961, "Relations nmong m sets of measures Among M Sets of Measure s " 
Psychometrika, 26, No.2, 129-149, 

- Caroll, J.D., 1968, "A generalisation Generalization of canonical correlation analysis 
Canonical Correlation Analysis to three -Three or more sets More Sets o f variables 
Variables ", Proc. 76 th Conv. Amer. Psych. Ass., 

- Kettenring, J.R., 1971, "Canonical analysis Analysis of several sets Several Sets o f 
variables Variables ", 58, 3, 433-450, or 

- Saporta, G., 1990, "Probability, analyse Analyse des donnees Bounces et 
statistique s Statistiques ", Technip, Paris. 

r0Q171 These methods however have certain limits. In particular, they do not 

allow to describe describing each set separately, the synthetic variables cannot be 
readily related to the various sets and, consequently, be physically interpreted. 
Furthermore, only the global inertia part they represent can be known, and not the 
inertia part they represent for each set. It is therefore not possible to filter each set 
independently. 

fQ0181 Another known analysis technique referred to as generalized principal 

components analysis (GPCA) allows to compare comparison of different sets of 
variables while describing each one. It is implemented on data of economic nature for 
example in : 

- Casin, Ph., 2001, "A generalization Generalization of principal component 
aralysfc Principal Component Analysis to K sets— Setg_o f variables Variables ", 
Computational Statistics & Data Analysis, 35, 417-428. 
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SUMMARY OF THE INVENTION 

fQ0191 The method according to the invention provides an approach wherein, by 

implementing a multivariable statistical technique applied to several sets of seismic 
records or logs available, obtained by exploration of an underground zone, the 
geophysical events readable on these sets of measurements are packed and filtered. The 
method comprises : 

- distributing these geophysical events into a certain number p of families of 
geophysical attributes having each a particular physical significances 

- forming, by combination of saM- the g eophysical attributes, synthetic variables in much 
fewer number than the number of sakUhe attributes, which are obtained by construction 
of an orthogonal vector base in each one of the p analysis sets consisting of the data of 
each one of said-the families, wherefrom is derived the formation of an orthonormal 
vector base describing the p analysis sets^; 

- using fekr-the orthonormal vector base for filtering and describing saM-the_geophysical 
events. 

[00201 The methodology for packing and filtering geo physical data in order to 

retrieve therefrom information on the nature of the sub soil is applied to packing and 
filtering sets of multi-domain geophysical data, i.e. either pre-stack seism ic data from a 
3D seismic survey grouped into a certain number of data cubes corresponding ea ch, for 
example, to the same offset (seismic trace cubes referred to as iso-offset cub es) or to the 
same angle of incidence (referred to as iso-incidence angle cubes), or seismic at tribute 
cubes obtained by seismic amplitude processing such as, f or example, elastic properties 
cubes from a joint stratigraphic inversion of pre-stack seismic data, these cubes being 
obtained bv implementing the method described for example in French paten t 2,800,473 
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(EP-1 .096.2701 filed bv the assignee. It then shapes the p r e-stack seismic information 



contained in seismic trace portions formed h v exploration of an underground zone, by 



means of new representative attribute s which are thereafter used for geologic 



interpretation or interpretation of an under g round reservoir. It can also apply to packing 
and filtering of pre-stack or post-stack 3 D seismic data from repeated surveys (4D 
seismic method) in order to interpret t h e fluid content variations of a petroleum 
reservoir under production. This method is a lso a pp li ca ble to packing of log , data of 
several physical natures by forming, fro m initial measurements, new log attributes in 
order to improve interpretation of the log data in ter ms o f lithologic and petrophvsical 
characteristics of the reservoir. 



r002*n F ormation of the synthetic variables is obtained for example by 

combining a principal components analysis technique or its derived techniques based on 
diagonalization of a norm (seeking eigenvectors ) and a simple regression technique. 
f00221 A ccording to an implementation mode, the synthetic variables can be 

subjected to a rotation so as to improve the correlations between initial variables, ife 
that is the geophysical attributes, and synthetic variables. 

r00231 T he synthetic variables can be formed from families comprising for 

example one or more seismic trace cubes obtained after a 3D seismic survey, each one 
corresponding to the same offset. 

r00241 Synthetic variables can also be formed from families consisting of 

seismic trace cubes obtained after a 3D seismic survey, each one corresponding to the 
same angle of incidence. 
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r00251 Synthetic variables can also be formed from families consisting of 

seismic attribute cubes formed by processing the seismic trace amplitudes obtained at 
the end of a 3D seismic survey. 

rpQ261 Synthetic variables can also be formed from families consisting of 

seismic trace cubes obtained after successive seismic surveys in the zone. 
r00271 It is also possible to form synthetic variables from families consisting of 

seismic attribute cubes formed by processing the seismic trace amplitudes obtained after 
successive seismic surveys in the zone. 

r00281 Synthetic variables can also be formed from families consisting of 

measurements of log attributes of different physical nature, such as the family of 
resistivity measurements, the family of measurements related to the formation 
radioactivity, the family of measurements related to the formation porosity, etc. 
rpQ291 The method is particularly advantageous in that it allows, when applied 

to seismic or log data : 

- to define new geophysical attributes, linear combinations of the initial geophysical 
attributes, summarizing the initial seismic or log data and allowing te-filtermg_of this 
data, while respecting the multidomain structure of the data ? : 

- to show if all the families of initial geophysical attributes, or only some of them, 
provide redundant information, and to measure the part of this redundant information^ 
and 

- to give a physical interpretation to the new geophysical attributes. 

r00301 The new geophysical attributes thus formed can then be analysed in order 

to retrieve the geologic properties of the seismic or log data. Two classes of methods 



can be conventionally used therefore : a first one using statistical estimation techniques 
and providing quantitative interpretation in terms of geologic properties, and a second 
one using pattern recognition techniques to analyse analyze t he seismic facies or the log 
facies (electrofacies) and provide qualitative information. 

[00311 Examples of these two types of methods are respectively implemented in 

the aforementioned pat e nts EP A- EP Patent 67 1,01 7_and (US[-] Patent 5,638,269) and 
FR-A -French Patent 2,768,818 and ft JS H Patent 6,051,651) filed by the 
applican t Assignee . 

BRIEF DESCRIPTION OF THE FK?URES DRA WINGS 

[00321 - Figures Figs, 1A to 1C show the organization of the seismic records in 

p sets : each set groups in a cube all of the trace portions corresponding to an identical 
physical criterion^: 

[0033] - Figures Figs, 2A to 2C respectively show the maps of the first three 

new seismic attributes Z (1) , Z (2) and Z (3) ? ; 

[0034] Figures Figs, 3A to 3D respectively show the first true temporal 

' plane (3 A) of the analysis window and the same plane reconstituted/filtered with 1 (3B), 
3 (3C) or 8 (3D) of the new seismic attributes for the cube of impedance P^i 

[00351 - Figures- Figs. 4A to 4D respectively show the first true temporal 

plane (4 A) of the analysis window and the same plane reconstituted/filtered with 1 (4B), 
3 (4C) or 8 (4D) of the new seismic attributes for the cube of impedance S ? ; 
[00361 Figures Figs. 5A to 5D respectively show the first true temporal 

plane (5 A) and the same plane reconstituted/filtered with 1 (5B), 3 (5C) or 8 (5D) of the 
new seismic attributes for the cube of density p?; 
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[00371 - Figures Figs. 6 respectively show the vertical sections, passing through 

well W2, true (6A), reconstituted/filtered with the first 8 new attributes (6B), and the 
differences (residues) (6C) between these two sections for the cube of impedance P ? : 

[00381 - Figures Figs. 7 A-C respectively show the vertical sections, passing 

through well W2, true (7 A), reconstituted/filtered with the first 8 new attributes (7B), 
and the differences (residues) (7C) between these two sections for the cube of 
impedance S ? ; 

[0039] - Figure Fig. 8 shows the location of the learning traces used to calibrate 

the classification function within the context of the supervised pattern recognition of the 
first 8 new seismic attributes^ 

[00401 - Figures -Figs. 9Ar and 9B respectively show the assignments of the 

seismic events and the most reliable assignments of these events within the context of 
the supervised pattern recognition of the first 8 new seismic attributes^ : and 

[00411 - Figures Figs. 10 At and 10B respectively show the assignments of the 

seismic events and the assignments of these seismic events within the context of the 
non-supervised pattern recognition of the first 8 new seismic attributes. 

DETAILED DESCRIPTION 

DETAILED DESCRIPTION OF THE INVENTION 

r00421 Several sets of geophysical records obtained by seismic exploration or 

logging in an underground zone studied are available. The geophysical events to be 
packed/filtered are characterized by geophysical attributes often organized in sets or 
families of particular physical significance (Figures 1A to 1C). These attributes can be 
the succession of the amplitude values of seismic trace portions or of elastic parameters 



12 



from a pre-stack joint elastic inversion, obtained for example by means of the method 
described in the aforementioned patent FR French Patent ,2,800,473 and_(EP[-l 
1,096,270) along the time window studied. These attributes can be, for a given depth, 
the measurements of a single set of logs obtained along a well, such as for example the 
values of the different resistivity tools, the values of the different tools measuring the 
porosity, etc. 

1 - Notations 

fQ043] Let X be a matrix consisting of p juxtaposed tables Xi,. . . Xj . . . X p : 

X = {X, ... Xj ... X p } , Xj = [Xy ... X jfk ... Xj, mj ] 
where Xj is a table consisting of mj variables X jjk of n individuals (in the statistical sense 
of the term). Matrix X groups together all the geophysical events available. These 
events are first organized in p sets or families Xj (Figure 1A). Each family Xj 
corresponds to a particular physical property. For analysis of the 3D seismic survey, this 
family has the form of a cube of values corresponding to this property, located by their 
spatial position (line, cdp, time/depth) (Figures IB and 1C). These cubes can for 
example consist of a set of amplitude values of n traces or trace portions associated with 
a common offset or a common angle or the spatial distribution of the values of a 
common attribute derived from these amplitudes such as, for example, the acoustic 
impedance in P or S waves, the density or Lame parameters, in a selected time window 
of size mj. In this example, each variable X j>k groups together the values of the seismic 
attribute considered (seismic amplitude or elastic property) for a temporal sample within 
the analysis window. In this example, a cube comprises as many individuals as traces 
available. For analysis of the logs, this family has the form of a sequence of mj series of 
values, located by their position along the well (time or depth). This sequence can 
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consist, for example, of all the logs of the tools of same physical nature, such as the 
various tools measuring the formation porosity, or the tools measuring various 
resistivity types, on a set of n time or depth points along one or more wells at the level 
of the underground formation studied. In this example, each variable X jjk groups 
together the values of the log attribute considered (resistivity, porosity, radioactivity,...) 
for a time or depth point within the analysis window. In this example, a sequence 
comprises as many individuals as time or depth points available along the well(s) 
analysed. 

[QQ441 The mj variables Xj, k of a table Xj define a space denoted by Sj and are 

assumed to be centred. Let D n be the diagonal matrix of the weights of the individuals. 

[0045] Let II j be an orthogonal projection operator on space Sj defined by : 

n j =x J (x<D n x J y l x<D tt 

2 - Calculation of the new attributes 

[00461 The purpose of Generalized Principal Components Analysis (GPCA) is 

both to describe the proximities between the p families of variables and to describe 
separately each set by means of synthetic variables that are linear combinations of the 
initial variables. Then, what is common to all the sets or families, or to some of them 
only, and consequently what is different can be shown. GPCA allows to reduce 
reduction of t he number of variables used to describe a data set and it allows 
multivariate filtering while keeping the multidomain structure of the initial data. To 
reach this objective, the first synthetic variable calculated, called Z (1) , has to be such 
that: 
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1. The p orthogonal projections Zj (1) of Z (1) on spaces Sj have to be as close as 
possible to one another; 

2. In each group j, the dispersion of the n individuals projected on Zj (l) is maximum. 
rOQ471 This is mathematically translated by the following conditions : 



1. max5> 2 (Z (, \z;'>) 

m J 

2. max£C0v 2 (Zj'\* y , v ) 

Z J v=l 

where R 2 is the multiple correlation coefficient and Zj (1) is the projection of Z (1) on 
space Sj normalized at 1 . 

[0048] A new synthetic variable Z (2) is then sought according to the same 

criteria. An orthogonality condition is added between the projections of variables Z of 
the same set. The procedure stops when an orthogonal base is found in each space Sj. 
The number of variables Z (k) does not exceed the largest dimension of spaces Sj. 

[00491 Rotation methods well known to the man skilled in the art can be 

implemented in each space thus found in order to optimize the correlation coefficients 
between the new variables and the initial variables. 

[00501 The first variable Z (1) is the normalized eigenvector of covariance matrix 

XX'Dn corresponding to the greatest eigenvalue. Z (1) is then projected on the p spaces Sj 
and the regression of the mj variables by Zj (1) is calculated in each space. The matrix of 
the p juxtaposed residues Xj (2) defines the new space studied. This new space 
corresponds to the initial space from which the influence of Zj (1) has been removed in 
each set Sj. The mj variables of table Xj (2) define a subspace Sj (2) of Sj. The second 
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component Z (2) is the normalized eigenvector of covariance matrix X (2) X (2)t D n of the 
residues corresponding to the greatest eigenvalue. Z (2) is projected on the p spaces Sj (2) . 
The regression of the mj variables by Zj (2) is calculated in each space Sj (2) . The 
regression carried out in each space Sj ensures that the projection of Z {2) on space Sj (2) is 
orthogonal to the projection of Z (1) on space Sj. 

fOQSn Calculation of the eigenvectors is carried out by means of any method 

known to the man skilled in the art. 

r00521 In practice, matrix XX l D n being of order n, where n can be very great, the 

normalized eigenvector V (k) of matrix X^DnX 00 associated with the greatest eigenvalue 
is calculated in a stage k. Vector Z (k) is then simply calculated by linear transformation 
Z°° = X^V 00 and normalized. 

[0053] Each set can be reconstituted with a small number of attributes. The 

initial variables of set j can be expressed as the sum of a linear combination of 
variables Zj w and of a residue. If all the variables Zj (k) of a set are used, the residue is 
zero. The coefficients of the linear combination are the coefficients calculated after each 
regression in spaces Sj. 

[0Q54] One of the properties of the method according to the invention is that 

variables Z 00 are orthogonal two by two. 

[DOSS] For interpretation of the synthetic variables calculated, we-then calculate 

for each set j: 

1. The square of the correlation between variable Z (k) and variable Zj (k) of set j, 
indicating the degree of proximity between these two variables; and 

2. the inertia percentage of set j explained by variable Zj (k) . 
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[00561 The correlations between variables Z (k) and initial variables X j>k are also 

calculated to show possible relations between the X j>k . The correlations between 
variables Zj (k) and Xj, k of the same set j, or the correlations between variables Zj (k) for 
the same k can also be calculated. 

fOQ571 The new synthetic variables thus created constitute the new geophysical 

attributes sought. 

3 - Application example 

fOQ581 The method according to the invention was applied to three cubes of 

elastic properties (impedance P denoted by Ip, impedance S denoted by Is and density 
denoted by p) obtained after pre-stack joint elastic inversion, by means of the method 
described in the aforementioned patent ^French Patent 2,800,473, of seismic records 
grouped in cubes of common incidence angles. Each elastic property generates a set of 
seismic attributes or variables. The analysis time window consists of 19 signal samples 
taken from seismic records. Wc consider the The_values of the elastic property for each 
temporal sample within the analysis window to be variables are considered. Each set 
thus consists of 19 variables made up of as many individuals as there are traces 
available. Since there are three variable sets, wc have in atotal 57 variables exists . The 
three parameter cubes have been flattened in relation to a reference surface. 
■ 00591 The synthetic variables Z 00 obtained summarize the respective 

information of the three elastic properties, while taking into account the relations 
between them. 
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3.1 - Analysis of the relations between impedances P, S and density p 

fQ0601 Table 1 hereunder gives the variance proportion in percentage explained 

by the projections of the first 8 synthetic variables Z (1) to Z (8) in the three sets. It shows 
that 8 variables are sufficient to summarize 87.4 % of the variance of impedance P, 94.4 
% of the variance of impedance S and 88.4 % of the variance of the density. 
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[00611 Table 2 hereafter gives the square of the correlation between synthetic 

variables Z (k) and their projections on the three sets. This value indicates the degree of 
proximity between the variable and the set. The closer the value is to 1 5 the closer the 
variable is to its projection and therefore to the set. 
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TOQ621 For example, Z (1) is very close to Is and close to p 5 but far from Ip, 

whereas Z (2) is close to Ip but very far from Is and p. Out of the first 8 variables Z (k) , 
three are close to Ip alone, two are close to both Is and p, one is close to Is alone, and 
two are close to p alone. It thus appears that Ip and Is do not contain the same 
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information and that part of the information contained in Is and in p seems to be 
common. 

3.2 - Filtering of impedances P, S and of the density 

r00631 The three elastic properties can be reconstituted with a limited number of 

new seismic attributes. The number of synthetic variables used for reconstitution can be 
different from one set to the next : it depends on the inertia percentage to be 
reconstituted/filtered. 

3.3 - Supervised and non-supervised pattern recognition of the new attributes 

[0064] It is possible to interpret geologically the new seismic attributes through 

a facies analysis with a supervised or non-supervised pattern recognition algorithm. 

1*00651 In supervised pattern recognition, it is necessary to select learning classes 

grouping together seismic events located in the vicinity of wells whose geologic 
interpretation has been carried out. This methodology allows to classify classification of 
all of the seismic events according to the geologic characters observed in the wells. 

[0066] In the present case, the geologic knowledge of the three available wells 

Wl, W2 and W3, and the general geologic knowledge of the zone studied lead us to 
considering three classes : a first class corresponding to wells Wl and W2 represents 
turbidite sands, a second class corresponding to well W3 represents turbidite muds and a 
third class represents hemipelagites. Wells Wl and W2 showing close geologic 
characteristics, only the seismic events located around well W2 are used, well Wl 
serving as a blind test. The location of well W3 being uncertain, the seismic events used 
to characterize class 2 are located around a virtual well defined in the proximity of W3. 
Finally, since we have t here are n o well for class 3, a virtual well is defined from our 
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general geologic knowledge. Each class consists of 121 seismic events (11 x 11). These 
events are used to calibrate a classification function which is then applied to the whole 
of the seismic events. 

rOQ671 Figures 8A. Figs. 8A and 8B respectively show the assignments obtained 

and the most probable assignments, i.e. having an assignment probability greater than or 
equal to 80 %. Well Wl, used in blind test, is correctly assigned to the same class as 
well W2. Well W3 is located at the boundary of classes 1 and 2. Class 2 is distributed 
along a channeling form. The map of Figure Fig. 8B shows that the results are stable : 
few events exhibit a bad assignment probability. 

[0068] In non-supervised pattern recognition, the learning classes are formed 

using the modes of a multivariate probability density function calculated from the whole 
of the seismic events, in a space generated by the first eight new attributes. The seismic 
events of higher probability density that constitute the peaks of this density function are 
selected to carry out learning of the classification function. This function is then applied 
to the whole of the seismic events considered to assign them to one of the classes : 
Figures 9 A. figs. 9A and 9B show the assignments obtained, as well as the most 
probable assignments, i.e. having an assignment probability greater than or equal to 
80 %. The maps obtained are very close to the maps of Figures 8A. Figs. 8A and 8B. 
Class 3 of the supervised analysis corresponds to classes 3, 4 and 5 of the non- 
supervised analysis. Wells Wl and W2 are assigned to the same class and well W3 is 
still located at the boundary between classes 3, and especially 1 and 2. This convergence 
of results validates the interpretation of the supervised approach. 



